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The energy-coupling factor (ECF) transporters are multi-subunit protein complexes that mediate uptake of 
transition-metal ions and vitamins in about 50% of the prokaryotes, including bacteria and archaea. Biological and 
structural studies have been focused on ECF transporters for vitamins, but the molecular mechanism by which ECF 
systems transport metal ions from the environment remains unknown. Here we report the first crystal structure of a 
NikM, 7YNikM2, the substrate-binding component (S component) of an ECF-type nickel transporter from Thermo- 
anaerobacter tengcongensis. In contrast to the structures of the vitamin-specific S proteins with six transmembrane 
segments (TSs), 7YNikM2 possesses an additional TS at its N-terminal region, resulting in an extracellular N-termi- 
nus. The highly conserved N-terminal loop inserts into the center of TfNikM2 and occludes a region corresponding 
to the substrate-binding sites of the vitamin-specific S components. Nickel binds to NikM via its coordination to four 
nitrogen atoms, which are derived from Metl, His2 and His67 residues. These nitrogen atoms form an approximately 
square-planar geometry, similar to that of the metal ion-binding sites in the amino-terminal Cu 2+ - and Ni 2+ -binding 
(ATCUN) motif. Replacements of residues in NikM contributing to nickel coordination compromised the Ni-trans- 
port activity. Furthermore, systematic quantum chemical investigation indicated that this geometry enables NikM 
to also selectively recognize Co 2+ . Indeed, the structure of 7YNikM2 containing a bound Co 2+ ion has almost no con- 
formational change compared to the structure that contains a nickel ion. Together, our data reveal an evolutionarily 
conserved mechanism underlying the metal selectivity of EcfS proteins, and provide insights into the ion-transloca- 
tion process mediated by ECF transporters. 
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Introduction 

The energy-coupling factor (ECF) transporters are 
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multi-subunit protein complexes that mediate uptake of 
transition-metal ions and vitamins across the membrane 
in many bacteria and archaea [1-7]. They consist of a 
substrate-specific integral membrane protein (S compo- 
nent, EcfS), a moderately conserved but ubiquitous trans- 
membrane protein (T component, EcfT), and homo- or 
heterooligomeric ABC ATPases (A components, EcfA). 
ECF transporters have been classified into two types. 
For type-1 transporters, all the components are encoded 
in the same operon, whereas for type-2 transporters, the 
EcfAA'T module is encoded in an operon, but the genes 
for S units are scattered around the genome. The S com- 
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ponents interact with the EcfAA'T module to form a ho- 
lotransporter complex. Previous structural studies have 
been focused on isolated S components, such as the ribo- 
flavin- [5], thiamine- [8], and biotin-specific [1] S units 
(RibU, ThiT, BioY), and on two type-2 holotransporters 
containing a folate-specific S component (FolT) [3] and 
an S component with predicted specificity for pyridoxine 
(PdxU, named HmpT in the original work [4]), respec- 
tively. All of these EcfS proteins have a similar structural 
fold with six transmembrane helices. The isolated S 
components were crystallized with their substrates that 
bind to a cavity facing the outside of the membrane. This 
observation correlates well with previous analyses of 
binding affinity that identified extremely low dissociation 
constants in the (sub)nanomolar range. 

ECF metal transporters belong to type-1 transporters. 
The mechanism by which ECF metal transporters selec- 
tively recognize metal ions remains largely unknown. As a 
unique feature of ECF metal transporters, function of the S 
component (NikM and CbiM for nickel and cobalt trans- 
porters, respectively) relies on additional small transmem- 
brane proteins represented by the "N" proteins in most 
cases or "K plus L" proteins for the rest. ECF-type nickel 
and cobalt transporters are the most common uptake sys- 
tems for these metals in prokaryotes [6, 9-11]. The vital 
importance of NikMNQO and CbiMNQO transporters (Q 
and O represent the T and A components, respectively) in 
the biology of microorganisms stems from the indispens- 
able roles of these two transition-metal ions in cells. 

Nickel is a cofactor of a variety of metalloenzymes 
involved in energy and nitrogen metabolism, carbon 
dioxide fixation and detoxification processes [9, 12-14]. 
Cobalt is a key constituent of vitamin B 12 [15]. Vitamin 
B 12 is the largest and the most structurally complicated 
water-soluble vitamin, and contains the biochemically 
rare element cobalt. In contrast to many archaea and 
bacteria that have the enzymes required for vitamin B 12 
synthesis, fungi, higher plants and animals are unable to 
synthesize vitamin B 12 . Therefore, microorganisms are 
the only sources of vitamin B 12 in nature [16]. Vitamin 
B 12 is involved in many critical biological processes in 
cells, such as DNA synthesis, fatty acid degradation and 
energy conservation [17, 18]. Vitamin B 12 deficiency 
can potentially cause severe and irreversible damage, 
especially to the brain and nervous system, and can also 
cause symptoms of mania and psychosis [19]. Although 
the biochemical mechanism for vitamin B 12 synthesis 
has been vastly studied [20, 21], the molecular mecha- 
nism by which the cells take up cobalt from the medium 
remains unclear. Several types of primary and second- 
ary active nickel and cobalt ion transporters have been 
reported [12], however in no case has the mechanism of 



selective metal binding and transport been elucidated. 

In the present study, we solved the crystal structure of 
7YNikM2, the S component of an ECF-type nickel trans- 
porter from Thermoanaerobacter tengcongensis . The 
overall structure of 7YNikM2 resembles the fold of the 
previously characterized S components of type-2 ECF 
transporters. It differs from them, however, by possess- 
ing an additional, N-terminal transmembrane segment 
that positions the N-terminus to the extracellular side of 
the membrane. The highly conserved N-terminal loop 
is crucial for the selection of metal ions. It inserts into a 
region corresponding to the substrate -binding sites of the 
vitamin-specific S components and contributes signifi- 
cantly to metal ion binding by donating three nitrogen at- 
oms. In total, four nitrogen atoms (the free amino group 
of Metl, the backbone nitrogen and the imidazole nitro- 
gen of His2, and the imidazole nitrogen of His67) form 
an approximately square-planar geometry for the specific 
recognition of Ni 2+ . Structural and quantum chemical 
analyses indicate that this architecture is also suitable for 
binding Co + , which is consistent with functional genom- 
ics analyses showing that the key residues for ion recog- 
nition are highly conserved among NikM/CbiM proteins 
from different species. Replacement of residues associ- 
ated with metal ion coordination reduces the Ni 2+ -uptake 
activity. Our structural studies, along with functional and 
theoretical analyses, reveal an evolutionarily conserved 
mechanism underlying the ion selectivity of this impor- 
tant class of transporters. 

Results and Discussion 

TtNikM2 is a nickel-specific S component 

Through functional genomics analyses, we identified 
two putative nickel-specific S components in T. tengcon- 
gensis: TTE0249 (7YNikM) and TTE1695 (7YNikM2) 
(Supplementary information, Figures S1-S3 and Table 
SI). Indeed, the X-ray fluorescence spectrum and the 
anomalous scattering signal screening results demon- 
strate that the metal ion present in the 7YNikM2 crystals 
is a nickel rather than a cobalt or iron ion (Figure 1A, 
Supplementary information, Figure S4 and Table S2). 
Using the nickel ion anomalous signal, we solved the 
structure of the 7YNikM2 protein with a C-terminal trun- 
cation (amino acids 1 to 209) at 1.83 A resolution (Sup- 
plementary information, Figures S4-S6). The full-length 
7YNikM2 (amino acids 1 to 230) structure was solved by 
molecular replacement method at 3.2 A resolution (Figure 
IB and Supplementary information, Table S3). The rmsd 
between the two structures is 0.25 A with 209 aligned 
Ca atoms. In the crystal, each asymmetric unit contains a 
single peptide and one nickel ion. The linker between the 
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Figure 1 Overall structure of TfNikM2. (A) The X-ray fluorescence spectrum as fitted shows multiple elements, e.g., CI, K, 
Ca, Fe, Ni and Zn ions. The nickel ion's Ka and Kp emission lines are shown as solid-dot purple lines. (B) Overall structure of 
TfNikM2. The TfNikM2 is shown in cartoon model, colored by a rainbow fashion. The Ni 2+ ion is shown as a sphere. The sec- 
ondary structural elements, N- and C-termini are indicated. (C) TfNikM2 is represented as an electrostatic surface model and 
two views related by 180° rotation along the vertical axis are shown. The blue, white, and red shadings represent positively 
charged, neutral, and negatively charged surface regions, respectively. 



a2 and a3 helices comprises two anti-parallel P-strands 
(Figure IB). The C-terminal end of 7YNikM2 forms a 
short a8 helix (amino acids 216-end). The surface of this 
a8 helix towards the cytoplasm is mainly hydrophobic, 
which suggests that this helix might interact with other 
components in the cytoplasm (Supplementary informa- 
tion, Figure S7). The overall structure of 7YNikM2 is 
hydrophobic with a positively charged belt facing the cy- 
toplasmic side of the membrane, which is in agreement 
with the positive-inside rule [22, 23] (Figure 1C). 

TtNikM2 has a unique N-terminus 

The overall structure of 7YNikM2 comprises seven trans- 
membrane segments with a short helix at the C-terminus 
(Figures IB and 2), and the additional transmembrane seg- 
ment at the N-terminus distinguishes it from the substrate- 



binding components of other ECF-type transporters [1-5, 
8]. The transmembrane segments 2-7 resemble the folds of 
RibU, ThiT and BioY, with rmsds of 3.27 A, 2.59 A, and 
2.88 A for 150, 133, and 140 aligned Ca atoms between 
r?NikM2 and the S units, respectively (Figure 2, left panel). 
The N-terminal loop (amino acids 1-9) in front of the al he- 
lix (amino acids 10-34) inserts into the center of 7YNikM2 
and occludes a region corresponding to the substrate-bind- 
ing sites of the vitamin-specific S components, ruling out 
the possibility that T?NikM2 recognizes any vitamins (Figure 
2, right panel). 

The unique N-terminus of NikM is required for substrate 
binding 

Previous studies have suggested that the amino group 
and the length of the N-terminal stretch of NikM/CbiM 
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proteins are essential determinants in metal recognition 
and transport [24]. The structure of 7YNikM2 reveals that 
the highly conserved first 9-residue loop inserts into the 
center of r?NikM2 and forms a series of hydrogen bonds 
with the surrounding residues (Figure 3B and Supple- 
mentary information, Figure S2). Interestingly, the free 
N-terminal amino group of Met 1 directly coordinates 
the nickel ion, and its side chain is stabilized by a deep, 
hydrophobic pocket that is mainly formed by Metll3, 
Alal50, Alal51 and Hisl90 residues (Figure 3B). The 
main chain amide nitrogen and the imidazole nitrogen 
of His2 also contribute to nickel binding. These three 
nitrogen atoms together with the side chain imidazole 
nitrogen of His67 form the nickel-binding site. This ob- 
servation is consistent with previous studies showing that 
adding an additional Met at the N-terminus or inserting 
an Ala between Metl and His2 resulted in the loss of 
transport activity [24]. Furthermore, the first 9 residues, 
and Gln94, GlylOl, Glyl02 and Asnl09 that interact 
with the N-terminal loop, are highly conserved in NikM/ 
CbiM proteins from different bacteria species (Figure 
3B, 3C and Supplementary information, Figure S2), sug- 
gesting a similar metal-recognizing mechanism among 
these proteins. 

Four nitrogen atoms form an approximately square- 
planar substrate-binding site 

The nickel ion is chelated by the nitrogen atom of the 




Figure 2 7fNikM2 has an additional transmembrane helix and 
an extra N-terminal loop compared to vitamin-specific S com- 
ponents. For TfNikM2, the extra N-terminal helix and loop are 
shown in blue and the rest of the protein is colored in cyan. 
RibU, ThiT, BioY are colored in wheat, slate, and lime, respec- 
tively. The right panel shows substrate-binding sites with 90° 
rotation to the left panel. The rioflavin, thiamine and biotin are 
colored in yellow, magenta and orange in sticks, respectively. 
The Ni 2+ is shown as a lightpink sphere. 



free N-terminal amino group of Metl, the main chain 
amide nitrogen and the N5 1 atom of His2, and the Ne2 
atom of His67 from the a2-a3 linker (Figure 4A). Un- 
like in other divalent metal transporters [11, 25, 26] and 
in nickel-binding proteins [12, 14], the four nitrogen 
atoms form an approximately square-planar geometry 
with slightly different distances from the nickel to the 
nitrogen atoms (Figure 4A and Supplementary informa- 
tion, Figure S8A). This arrangement is reminiscent of 
the approximately square -planar geometry of the amino- 
terminal Cu 2+ - and Ni 2+ -binding (ATCUN) motif, which 
has been found in certain species of albumins (e.g., HSA, 
BSA, and RSA), human sperm protamine P2a, neurome- 
dins C and K, and histatins [27-35]. Unlike the metal ion 
binding mode in 7YNikM2, the ATCUN motif consists of 
the N-terminal amine, the backbone amide nitrogens of 
the second and third residues, and an imidazole nitrogen 
of His4 (Supplementary information, Figure S8). Consis- 
tent with our structural observation, the nickel AT-edge X- 
ray absorption near-edge structure (XANES) results also 
indicate a square-planar geometry around the nickel ion 
of 7YNikM2 in solution. The shoulder feature at 8339.5 
eV present in 7YNikM2 and in Ni(dmg) 2 is attributed to 
the ls-to-4p z transition [36], indicative of a characteris- 
tic of the almost square-planar coordination geometry 
around the nickel (Figure 4B). 

The nickel-chelating residues are stabilized by a series 
of hydrogen bonds (Figure 4A). The hydrogen atoms 
from the N-terminal free amine form hydrogen bonds 
with the carbonyl group of Thr65 and a water molecule 
(Figure 4C). This water molecule forms additional hy- 
drogen bonds with the carbonyl group of His67 and the 
Ne2 atom of Hisl90. The N51 atom of Hisl90 forms a 
hydrogen bond with the carbonyl group of Metl 86. The 
Ne2 atom of His2 forms a hydrogen bond with the side 
chain of Gln94. The N51 atom of His67 forms a hydro- 
gen bond with the carbonyl group of Ala68, and the side 
chain of Glul98 forms another hydrogen bond with the 
backbone amide of Ile69 (Figure 4C). Interestingly, the 
Ni 2+ has two positive charges, but from the structural 
analyses we predict that only one hydrogen from the 
main-chain amide of His2 could be deprotonated. We 
propose that this elaborate hydrogen bond network not 
only stabilizes the planar nickel ion-binding geometry, 
but also contributes to the charge neutralization. More 
interestingly, these residues involved in the formation of 
the hydrogen bond network are highly conserved among 
the NikM/CbiM proteins from different species (Supple- 
mentary information, Figure S2), suggesting that the 
mechanism of nickel or cobalt recognition is most likely 
conserved. 

To confirm the biological relevance of the structural 
observations, we attempted to examine Ni 2+ uptake in 
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Figure 3 The N-terminus of 7YNikM2 is responsible for Ni 2 * binding. (A) The N-terminal loop of 7fNikM2 inserts into a deep 
pocket. The a1 -helix and the N-terminal loop are shown in ribbon representation and the rest of the protein is shown in 
an electrostatic potential surface. The blue, white, and red shadings represent positively charged, neutral, and negatively 
charged surface regions, respectively. The Ni 2+ ion is shown as a sphere. (B) A close view of the N-terminal substrate-binding 
pocket. The residues from N-terminal loop and from the rest of the protein are shown in yellow and magenta sticks, respec- 
tively. The water molecule is shown as a red sphere. The Ni 2+ -nitrogen coordinate bonds and hydrogen bonds are shown as 
green and red dashed lines, respectively. (C) The sequence alignment of the highly conserved N-terminus (amino acids 1-9) 
of CbiMs and NikMs from different species. 



recombinant Eischerichia coli cells producing 7YNik- 
M2N2. It was previously shown that M plus N proteins 
(in the absence of Q and O components) are the minimal 
unit of ECF-type metal transporters required to medi- 
ate basal levels of substrate uptake [6, 24]. However, 
we failed to detect any metal ion uptake above back- 
ground. It is possible that J?NikM2 and 7YNikN2 from 
the thermophilic T. tengcongensis do not assemble into 
a functional transporter or it is likely that other compo- 
nents are required for the establishment of a functional 
transporter in the mesophilic E. coli host. Moreover, M 
and N proteins were previously shown to only loosely 
interact with each other [24]. Therefore, we used the 
homologous nickel transporter, 7?cNik(MN)QO, from 
the purple bacterium Rhodobacter capsulatus to test the 
biological relevance of the structural observations (Fig- 
ure 4D). In this system, the NikM and NikN proteins are 
fused naturally to form a single Nik(MN) peptide. Given 
the sequence similarity between the nickel transporters 
from T. tengcongensis and R. capsulatus (Supplemen- 



tary information, Figure S9), we replaced the structur- 
ally important residues in the NikM part of i?cNik(MN) 
QO and tested the effects on nickel-uptake activity. 
Replacement of His67 (i?cNik(M His67Ala N)QO), which is 
involved in nickel binding, abolished transport activity. 
Likewise, replacements of residues stabilizing the planar 
geometry greatly (Hisl94Ala, Hisl94Glu, Glu202Gln, 
Glu202Asn, Glu202A) or at least significantly (Thr65Ala) 
impaired the Ni~ + -uptake activity, whereas replacements 
that preserve hydrogen bonds (Gln94Glu, Hisl94Gln, 
Glu202Asp) exhibited only a moderate effect (Figure 
4D). The latter finding is consistent with our prediction 
that the hydrogen bond network contributes to the charge 
neutralization (Figure 4D). The essential roles of Metl 
and His2 for the Ni 2+ -uptake activity have been reported 
[6, 24]. 

The substrate-binding site dictates the ion specificity 

Structural analysis suggests that only the hydrogen 
atom from the backbone amide of His2 can be deproton- 
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Figure 4 The approximately square-planar geometry of Ni 2 *-binding site. (A) The detailed binding mode of the Ni 2 *. The resi- 
dues involved in Ni 2+ coordination and stabilization of the binding pocket are shown as sticks. The Ni 2+ -nitrogen coordinate 
bonds and hydrogen bonds are shown as green and red dashed lines, respectively. (B) Comparison of the Ni K-edge XANES 
spectra of the NiO, Ni(dmg) 2 in solid and TfNikM2 in solution. The observed distinct shoulder feature at 8339.5 eV present in 
the TfNikM2 and Ni(dmg) 2 is attributed to the 1s-to-4p z electronic transition in a square-planar metal complex. The structures 
of the Ni 2+ -binding sites of NiO, Ni(dmg) 2 and TfNikM2 are shown in sticks. (C) The electron distribution of the Ni 2+ -binding 
pocket. The residues are shown in stick model, with hydrogen atoms displayed. The distributions of electrons are indicated 
by 8" and 8 + . (D) Ni 2+ -uptake activity of E. coli cells containing the RcNik(MN)QO transporter from R. capsulatus with wild-type 
or mutant RcNikM. Each assay was performed in triplicate. Mean values ± standard deviations are shown. Numbers in paren- 
theses correspond to residue positions in RcNikM, where the numbering differs between RcNikM and HNikM. 



ated upon coordination of Ni' + in the cavity formed by 
four N-atoms (Figure 5A). Ni 2+ has two positive charges 
while the deprotonated cavity possesses one negative 
charge. It is thus imperative to understand how the pro- 
tein neutralizes the extra positive charge and how this 
metal-binding geometry dictates ion selectivity. To eluci- 
date the nature of the selectivity of the 7YNikM2 protein 
for divalent nickel ion, we have performed systematic 
quantum chemical investigations of the chemical bond- 
ing, charge distribution and binding strength of the pro- 
tein with a series of biologically relevant divalent metal 
cations (M = Ca 2+ , Mg 2+ , Fe 2+ , Co 2+ , Ni 2+ , Cu 2+ and Zn 2+ ) 
(Figure 5B and 5C). We selected a fragment L (Figure 



5A) of the local structure of the protein and the dangling 
bonds are saturated with hydrogen atoms to fulfill the 
correct valence states of each atom. As the protein is de- 
protonated upon coordination of the metal ions, the bind- 
ing ligand is negatively charged. The Gibbs free energies 
for M 2+ + L~ — > ML + were calculated as AG = G(ML + ) - 
G(M 2 ^) - G(L~), where the entropy effects are included 
along with the quantum mechanical-binding energies. 
The negative energy indicates the binding is thermody- 
namically exothermic. The metal-binding constants K p 
were calculated from AG = — RT lnAT p . The energies of the 
M 2+ and ML + ions are determined with all possible multi- 
plets and spin-multiplicities. The lowest-energy electron 
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Figure 5 The approximately square-planar geometry of the metal-binding site dictates substrate specificity of 7fNikM2. (A) 
Structure of the 7fNikM2 nickel coordinate cavity. The nitrogen atoms are defined as N1, N2, N3, and N4. (B) Binding energy 
of M 2+ + L~— > ML + (M represents different divalent metal ions) as calculated with BHaH functional. (C) Splitting of the 3d or- 
bital in NiL* and 3D contours (isovalue = 0.03) of Kohn-Sham molecular orbitals. The calculated Kohn-Sham energy levels 
and the contour diagrams of the frontier molecular orbitals of NiL + are shown, where Ni 2 * has a d 8 electron configuration. (D) 
Structural alignment of the TfNikM2-nickel and 7fNikM2-cobalt structures. The structures with nickel and cobalt are shown in 
magenta and slate, respectively. The residues binding the metal ions are shown in sticks. 



configurations are used as the references for the Gibbs 
free energies. The atomic electron configurations of M 2+ 
used are consistent with those proposed by Johnson et al. 
[37]. 

The density functional theory (DFT) calculation [37- 
39] results indicate that when the divalent metal cations 
bind to the substrate -binding cavity, Ca + and Mg + have 
the least binding energy due to the pure ionic interaction 
(Figure 5B, and Supplementary information, Table S4), 
while the transition-metal ions have larger binding ener- 
gies due to both ionic and covalent interactions between 
the N-atoms and the metal 3d-orbitals (see the contours 
of the NiL + orbitals in Figure 5C). In the structures, li- 
gand field theory (LFT) predicts that the d-orbital mani- 
fold of the metal ions will split into five groups, with one 



strong anti-bonding orbital (d X 2-y2) lying above the other 
four d-orbitals (Figure 5C). This ligand field splitting 
pattern indicates that a d 8 configuration should have larg- 
er binding energy than the others. Indeed, the DFT cal- 
culations support this qualitative conclusion, with Co 2+ 
(d 7 ) and Ni 2+ (d 8 ) having Gibbs free energies of -22.47 
eV and -23.09 eV (electron volt), respectively (Figure 
5B and Supplementary information, Table S4). The Cu 2+ 
(d 9 ) ion tends to have even larger free energy, but this ion 
might prefer a different local structure due to the Jahn- 
Teller effect. 

The calculated atomic polar tensor (APT)-derived 
charges [40] indicate that the four N-atoms have dif- 
ferent charges (Supplementary information, Table S5). 
Their charges are not only consistent with the distance 



www.cell-research.com | Cell Research 



@ Crystal structure of NikM 



difference between the Ni-N„, (m = 1, 2, 3, 4) coordina- 
tion bonds (Figure 5A and Supplementary infonnation, 
Figure S8A), but also consistent with the notion that the 
surrounding hydrogen-bond network contributes to the 
charge neutralization (Figure 4C). These data suggest that 
the metal ion type and the metal-ligand bonding strength 
might constitute the most important elements to deter- 
mine the metal ion selectivity. To confirm the theoretical 
results, we purified and co-crystallized r?NikM2 1 _ 209 with 
CoCl 2 and solved the structure at 2.5 A resolution (Sup- 
plementary information, Figure S10 and Table S3). The 
structure has a 0.22 A rmsd with the structure that con- 
tains nickel ion and has almost no conformational chang- 
es (Figure 5D), which suggests that 7YNikM2 could bind 
both nickel and/or cobalt ions. These results confirm the 
hypothesis that CbiM and NikM proteins utilize similar 
mechanisms to recognize their substrates. They are also 
consistent with the previous biological observation that 
the CbiMNQO complexes of both Salmonella enterica 
serovar Typhimurium and R. capsulatus not only import 
cobalt ions but also transport nickel ions [6]. 

Conclusions 

ECF transporters have evolved into two types with 
high selectivity for either vitamins or transition-metal 
ions. Our studies unravel that the specific Co 2+ and/or 
Ni 2+ recognition by NikM is achieved by the presence of 
an additional N-terminal transmembrane helix that en- 
ables the very N-terminal end of the protein with its three 
metal ligands to deeply insert into the substrate-binding 
site, thereby preventing entrance of larger molecules and 
establishing ion selectivity. This represents another fasci- 
nating feature of the structurally and mechanistically un- 
usual ECF transporters. The mechanistic studies of how 
these protein complexes exert their micronutrient-trans- 
port functions could be valuable for the development of 
new antibiotics since ECF transporters are essential for 
pathogenic bacteria. 

Materials and Methods 

Genome-wide analysis of phylogenetic relationships of 
Nickel and Cobalt transporters 

The cbiM and nikM genes were obtained from KEGG (the Kyo- 
to Encyclopedia of Genes and Genomes) ORTHOLOGY Database 
(KEGG ORTHOLOGY Entry K02007) (Supplementary informa- 
tion, Table SI). The collected sequences were aligned using the 
MUSCLE (multiple sequence comparison by log-expectation) pro- 
gram [41] and further adjusted manually using BIOEDIT. The op- 
timal substitution model of amino acid substitution (RtREV+G+F) 
was selected using the modelGenerator version 0.84 program 
[42]. The phylogenetic tree in Supplementary information, Figure 
SI was constructed using the maximum likelihood method with 
PhyML software [43]. One hundred bootstrap replicates were per- 



formed to obtain the confidence support. 

The X-ray fluorescence spectrum assay 

The microscopic X-ray fluorescence (u-XRF) measurements 
are performed at Beijing Synchrotron Radiation Facility (BSRF) 
4W1B endstation, which runs 2.5 GeV with current from 150 
mA to 250 mA. The X-ray beam impinged on 7/YNikM2 crystal 
is monochromatized at 15 keV and is focused down to 50 um in 
diameter by the poly-capillary lens. The detection limit of the ex- 
perimental setup is ~ Hg/g. The 7HSfikM2 crystal is pasted on XRF 
tape (TF-500) at room temperature and then attached to a frame- 
shape sample holder. The Si (Li) solid state detector is used to de- 
tect X-ray fluorescence emission lines with live time of 100 s per 
point. Several XRF spectra were collected at the sample and the 
blank tape, respectively. The data reduction and process were per- 
formed using PyMCA package with linear polynomial background 
accounting for scattering [44]. Because the X-ray might penetrate 
the sample and strike the tape, the XRF spectrum of pure 7YNikM2 
crystal is achieved by subtracting the XRF spectrum of blank tape. 

Protein preparation 

PCR-amplified full-length nikM2 from T. tengcongensis was 
subcloned into pET21b vector with a TEV protease site before the 
His-tag. Overexpression of TtnikM2 was induced in E. coli BL21 
(DE3) by 0.5 mM isopropyl-[3-D-thiogalactoside (IPTG) when the 
cell density reached OD 600 = 1.0. After growth in Terrific Broth for 
3 h at 37 °C, the cells were harvested, homogenized in buffer con- 
taining 20 mM Tris-HCl (pH 8.0), 500 mM NaCl, 20 mM imid- 
azole, 1 mM PMSF and disrupted by sonication. Cell debris were 
removed by low-speed centrifugation for 15 min. The supernatant 
was collected and applied to ultracentrifugation at 150 000* g for 
1 h. The membrane fraction was harvested and suspended in 20 
mM Tris-HCl (pH 8.0), 500 mM NaCl, 20 mM imidazole, 1 mM 
PMSF. After the addition of Triton X-100 at a final concentration 
of 1% (v/v), the sample was incubated for 10 h with slow stirring 
at 4 °C. After another ultracentrifugation step at 150 000* g for 1 h, 
the supernatant was collected and loaded onto Ni 2+ -nitrilotriacetate 
affinity resin (Ni-NTA) and rinsed with buffer containing 20 mM 
Tris-HCl (pH 8.0), 500 mM NaCl, 40 mM imidazole and 0.02 % (w/ 
v) DDM. The protein was eluted from the affinity resin with buf- 
fer containing 20 mM Tris-HCl, pH 8.0, 500 mM NaCl, 300 mM 
imidazole and 0.02 % (w/v) DDM. The protein was concentrated 
to 5 mg/ml before further purification by gel filtration (Super- 
dex-200 column, GE Healthcare) in buffer containing 20 mM Tris- 
HCl (pH 8.0), 500 mM NaCl, 2 mM DTT, 0.02% (w/v) DDM. 
The peak fractions were collected and concentrated to 3.2 mg/ 
ml. 7YNikM2i_ 209 (amino acids from 1 to 209) was cloned, overex- 
pressed and purified in the same way. The 7YNikM2 1 _ 2 o9 with the 
divalent cobalt ion was purified by same method but using the Co- 
NTA column instead of the Ni-NTA. 

Ciystallization, data collection and structural determination 
Extensive crystallization trials were performed for the purified 
7/YNikM2 proteins. Crystals were grown at 18 °C by the hanging- 
drop vapour diffusion method. After optimization, the full-length 
7YNikM2 crystals were grown by mixing 1 (4.1 of protein (3.2 mg/ 
ml) with 1 ul of reservoir solution (38% w/v PEG300, 0.05 M 
glycine, pH9.5, 0.1 M NaCl, 0.2 (ri 35% w/v 1,6-hexanediol, 
0.2(4.1 6.6% w/v cyclohexylethanoyl-N-hydroxyethylglucamide). 
77NikM2 crystals appeared and grew to full size in about one 
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week. The 7 , iNikM2 I _ 20 9 crystals were grown in the same way ex- 
cept that the concentration of PEG300 in the reservoir solution is 
35.5% w/v. The r/NikM2,_ 109 crystals with cobalt ion were grown 
in the same way except that the buffer changed to MES with the 
pH 6.4. All the crystals were frozen in liquid nitrogen directly to 
collect the data. 

All the crystal data were collected at the Shanghai Synchrotron 
Radiation Facility (SSRF) BL17U, integrated and scaled using 
the HKL2000 package [45]. Further processing was carried out 
using programs from the CCP4 suite [46]. Data collection sta- 
tistics are summarized in Supplementary information, Table S2 
and Table S3. The nickel sites were located using SHELXD [47] 
from Bijvoet differences in the Ni-SAD data. Heavy atom posi- 
tions were defined and the phases were calculated with PHASER'S 
SAD (single anomalous diffraction) experimental phasing module 
[48]. The real-space constraints were applied to the electron den- 
sity map in DM. The final model rebuilding was perfonned using 
COOT [49] and the protein structure was refined with PHENIX 
[50] using NCS and stereochemistry information as restraints. The 
structure of full-length 7flSfikM2 was solved by molecular replace- 
ment method with the 7YNikM2 ,_,(><, structure as the search model. 
All structural figures were generated in PyMOL [51]. 

XANES data collection 

The concentration of the 7YNikM2 was about 0.5 mM. Samples 
were flash frozen in liquid nitrogen. During the XANES measure- 
ment, the sample was maintained at a temperature of 80 K. All 
measurements of the Ni AT-edge XANES were carried out on beam 
line BL14W1 at the Shanghai Synchrotron Radiation Facility 
(SSRF, China) in the fluorescence yield mode using a double-crys- 
tal Si(l 1 1) monochromator and a 32-element array of Ge detectors. 
A N 2 -filled ionization chamber was used to measure the incident 
flux. The storage ring was working at 3.5 GeV and experiments 
were performed with an electron current 220 mA. Simultaneous 
energy calibration was performed by measuring the absorption of 
Ni foil. 

The XANES spectrum is highly sensitive to the coordination 
geometry around the absorbing atom. A shoulder peak appeared a 
few eV below the main absorption edge in the XANES spectrum 
which is attributed to the 1 s-to-4p z electronic transition is a signa- 
ture of a square-planar complex [36]. The presence or absence of 
the ls-to-4p z transition has been used to predict the coordination 
geometry in several model compounds [52, 53]. 

Construction of RcNik(MN)QO transporter variants and 
transport assay 

Plasmid pRcNik(MN)QO encoding the natural Nik(MN) fusion, 
NikQ and a C-terminally FLAG-tagged NikO [6] was the basis for 
all constructions. Individual mutations were introduced by two 
rounds of PCR. In the first round, amplicons were generated with 
a mutagenic forward or reverse primer and a corresponding primer 
directed against the 3 '-end of nikMN or the upstream region in the 
vector, respectively. Those amplicons were used as primers in a 
second round of PCR together with the appropriate non-mutagenic 
primer of the first round. The second-round products were treated 
with appropriate restriction endonucleases and were used to re- 
place the wild-type fragments. All synthetic DNAs were verified 
by nucleotide sequencing. Mutant plasmids were introduced into 
E. coli XL 1 -Blue and nickel-uptake activity was determined as 
described [6]. Briefly, cells were grown in lysogeny broth (LB) 



over night, diluted 1: 100 in LB supplemented with ampicillin (100 
(ig/ml), IPTG (1 mM) and 63 NiCl 2 (500 nM), and incubated under 
shaking at 37 °C for seven hours. Cells were harvested, washed 
twice with Tris-hydrochloride (50 mM, pH 7.5) and concentrated. 
Metal accumulation was quantitated by liquid scintillation count- 
ing and is expressed as pmol/mg protein. 

Quantum chemical analysis 

The theoretical calculations were performed with Gaussian 09 
(G09) program [54]. We used the LANL2DZ effective core poten- 
tials and corresponding valence basis sets for the metals [55] and 
6-311+G* all-electron basis sets for C, N, O, and H [56]. Because 
of the complexity of 3d-transition metals, we investigated the 
binding energies using a variety of approximate exchange-corre- 
lation (XC) functionals as benchmark. The XC functionals used 
include the generalized gradient approximation (GGA) of Perdew- 
Burke-Ernzerhof (PBE) [57], hybrid GGA of PBE0, B3LYP, and 
BHaH (Becke's half-and-half hybrid functional), highly nonlocal 
hybrid meta-GGA functional M06-2X and full Hartre-Fock func- 
tional M06-HF [58]. The basic trends of the calculated energies 
and charges are similar from these XC functional. Based on the 
charge and spin-density distribution, the performance of the BHaH 
functional is satisfactory and we only present results calculated 
from this functional. The results are also compared with those cal- 
culated using Amsterdam Density Fucntional (ADF) program [59]. 

The vibrational frequencies were calculated by using harmonic 
approximation. The zero-point energy (ZPE) is included in the 
calculated binding energies. The atomic polar tensor (APT) de- 
rived charges were calculated to provide the information on charge 
distribution upon different ligands binding [40]. The energies cal- 
culated are listed in Supplementary information, Table S4 and the 
atomic net charges are listed in Supplementary information, Table 
S5. 

PDB deposition 

The coordinates and diffraction data of 77NikM2i_ 209 (with 
Co 2+ ), 7YNikM2,_ 209 (with Ni 2+ ) and the full-length 7YNikM2 with 
Ni 2+ have been deposited into the Protein Data Bank with acces- 
sion code 4M5C, 4M5B and 4M58, respectively. 
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